LETTER TO THE EDITOR 



The noncommutative degenerate electron gas 



F S Bemficaf and H O Girottif 

f Instituto de Ffsica, Universidade Federal do Rio Grande do Sul, Caixa Postal 
15051, 91501-970 - Porto Alegre, RS, Brazil 

E-mail: f bemf icaOif .uf rgs .br, hgirottiOif .uf rgs .br 

Abstract. The quantum dynamics of nonrelativistic single particle systems involving 
noncommutative coordinates, usually referred to as noncommutative quantum 
mechanics, has lately been the object of several investigations. In this note we 
pursue these studies for the case of multi-particle systems. We use as a prototype 
the degenerate electron gas whose dynamics is well known in the commutative limit. 
Our central aim here is to understand qualitatively, rather than quantitatively, the 
main modifications induced by the presence of noncommutative coordinates. We 
shall first see that the noncommutativity modifies the exchange correlation energy 
while preserving the electric neutrality of the model. By employing time-independent 
perturbation theory together with the Seiberg-Witten map we show, afterwards, that 
the ionization potential is modified by the noncommutativity. It also turns out 
that the noncommutative parameter acts as a reference temperature. Hence, the 
noncommutativity lifts the degeneracy of the zero temperature electron gas. 
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The first paper on quantum field theories formulated in a noncommutative space- 
time manifold was published in 1947 |TJ, although the idea that a noncommutative 
space-time manifold might provide a solution for the problem of ultraviolet divergences 
seems to have been suggested long before jSJ. The subject was, perhaps, abandoned due 
to the success of renormalization theory and its revival is rather recent and related to 
string theory. Indeed, the noncommutative Yang-Mills theory arises as a limit of string 
theory [3| and it was extracted by Seiberg and Witten jl] by starting from the open 
string in the presence of a magnetic field. More details on this and related subjects 
can be found in the already existing review articles El El IH1 13 UD] and also in the 
specialized literature. 

On the other hand, noncommutative quantum mechanics has also been under 
scrutiny [HI [TJl HHj . The main outcome, in the case of single particle systems, is that 
a modification of the equal-time algebra obeyed by the basic position observables acts 
as a source of new interactions which may or may not preserve the original symmetries. 
The present paper is dedicated to study the physical consequences of introducing 
noncommutative coordinates in the case of quantum many-particle systems. 

We consider, as a prototype, an idealized high density degenerate electron gas 
occupying a volume V = L 3 . As it is currently assumed, the electrons (electric charge 
— e) are in the presence of a uniform background of positive ions (electric charge +e) 
that makes the whole system electrically neutral. For this to be the case the number 
of electrons must equal the number of ions (N). The ions are much heavier than 
the electrons and will, then, be considered as static. Although this system has been 
extensively described in textbooks we make a small digression here to pinpoint its 
highlights. By election, the degrees of freedom of each electron (E) are the Cartesian 
positions {X J a E } and linear momenta {P ] a E }, together with the spins {S 3 a E }. For the 
ions (B), the corresponding observables are, respectively, {X J a B }, {P 3 a B } and {S J a B }. 
Small case letters from the beginning of the Latin alphabet (a, ft,...) designate the 
particle, while small case letters from the middle of the Latin alphabet . . .) only 
run from 1 to 3 and identify the Cartesian component of the corresponding observable. 
Observables associated with the electrons commute with those associated with the ions. 
The phase space electron degrees of freedom obey the standard equal-time algebra 

[K,b,H,e] = MPS*, (1) 
[Pa , e > Pi, e] = • 

We emphasize that all position observables commute among themselves. The equal-time 
algebra for the ion phase space variables can be obtained from Eq.((TJ) just by replacing 
E by B. The algebra of the spin components will not be explicitly displayed. 
Structurally, the more general form of the total Hamiltonian reads 



H = He + H B + H EB 



(2) 
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The Hamiltonian H B describes the free dynamics of the electrons plus the Coulomb 
interaction among them. Hence, in the position representation (X 1 \x) = 
x l \x) , P l — > p' 1 = — ihd/dx l ), 

Nii -, N 

a=l a^b 

where 

V(|f|) = e 2 ^ (4) 
\r\ 

and /i is a damping factor needed to secure, in the thermodynamic limit, the convergence 
of each term in the right hand side of Eq. (J2J) 14 . As for H B one writes 



1 N 

H B = 77 ^ V (l ^a,B - ^6,B |) 



V 

2 

^ y d 3 x„ jjB J d 3 x b:B n(x at B)n(x b:B )y(\xa,B - x bfB \)- (5) 

The absence of kinetic term in Eq.fjHJ) reflects the fact that the ions are static. 
Furthermore, the continuous nature of the ion background is taken into account by 
replacing the discrete summations by continuous integrals. This brings into play the 
new variable, n(x a>B ), known as the ion density. Finally, H BB is taken to be 

N N 



H E B = ~ V(| Xg,E - %b,B [) 

a=l 6=1 
N 

-» / ^ 3 x bjB n(x btB ) V(\x a>B - x b>B \). (6) 

o=l J 

We shall always be working in the approximation n(x a ^ B ) = constant = N/V. 
It has long been shown [T3] that H, in Eq.©, can be cast 



H = H + Hj , (7) 



where 



Eft k + / \ 

"2^T C L C £A' (8a) 

fcA 

F ' = Y E'E^W.A.^A.^A,^^. (86) 
kpq A l A 2 

Here, ct (c^ A ) are the creation (annihilation) operators of electrons of momentum k 
and spin A, whereas q = \q\. Furthermore, the prime in the summation symbol implies 
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that the momentum q = is excluded. Within the framework of time-independent 
perturbation theory, the main outcomes, including contributions up to the second order, 
may be summarized as follows. The unperturbed ground state (Fermi) energy (E^) is 
given by [13] 

„2 



E, 



(0) 



N 



2.21 



(9) 



2a r% 

where a is the Bohr radius, r s = r /a Q and |vrrQ = V/N. Moreover, its first and second 
order perturbative corrections were, respectively, found to read [131 E3 EH H7j 

„2 



E, 



(i) 



2a n 



N 



0.916 



(10) 



and 



E 



(2) 



2a n 



N 



,(2)r 



.(2)6 



AT [0.0622 lnr s - 0.094] . 

2ao 



Here [pS], 



e ( 2)6 = Af^f d s k / d 3 p e(l-k)e(l- P ) 

167r>,/ q 2 7| £+g | >1 ./ M>1 [q-(q + k _ p)](q + k _ p)2 

= 1^2- C(3) « 0.048, (12) 

is the exchange correlation energy in Rydberg's units. We shall designate by £(x) the 
Heaviside step function. In order to work with dimensionless vectors (p) we define 
kp~p = p, where k E = r^ 1 is the Fermi wavenumber. Also, p = |p|. Here ends 

our brief summary about the degenerate electron gas. 

We turn next into studying the implications of replacing X l a E — > Q l a E , P % a E — > 
P % a E , now obeying the equal-time phase space algebra 



[Qa , E > Qb , e] - 2 i 5a , b > 

[Ql E ,PiJ = ih^S ab , (13) 

[P l a, E ,Pi >E ] = 0. 

The distinctive feature of the new position observables (Q l a E ) is that they do not 
commute among themselves. This lack of noncommutativity is characterized by the real 
antisymmetric constant matrix (0%). An explicit representation for this algebra has 
already been obtained [TTJ [T2 [THl UU| after realizing that (see Eqs.([TJ) and JH1) 

Qi, E = k >e - i®% p i,E- (14) 

A similar modification should be introduced for the ions. However, for static ions Eq.(|14j) 
reduces to Q\ B = X l a B . 
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As for the Hamiltonians, the replacement X l a E 
H — > TC, such that 



Qa,Ei Pa,E 



where 



T~t — 7~(-e + 7~t b + 7~Leb 



P l a E amounts to 



(15) 



Hi 



N i i 

V I'a.E 



a=l 



2m 



1 - 



'>a,E ~ <Pb,E 



(16) 



and 



T-Cb = ^ J d 3 x a , B J d 3 x biB n(x a)B )n(x b)B )V(\x ajB - x b:B \) 



H 



EB 



N „ 

^ / d 3 x b , B n(x bjB ) V ( 

a=l J 



^>a,E — X b B 



For simplifying purposes, we have introduced the notation 



, E — X a , E ^ ® J Pa , E • 



(17) 



(19) 



As already stated, {x l a E } denotes the set of eigenvalues of the operator X % a E , whereas 
p l a E = —ihd/dx l a E represents P % a E in the basis defined by the common eigenvectors 
Of 



We next focus on V 



*>a,E - %b,B 



when acting on an arbitrary but differentiable 



function ^{x l a E , x\ B ). By taking into account Eq.([19|) one finds that 



V 



x 



a,E 



- ^ Q ^Pa,E ~ X 



i 

b,B 



*K,e,4,b) 



(2 



-\— f d 3 A;V(fc)e iA;l{ <^-i e ^^^-^.s)^(x J a 

T^Y' J 

d 3 kV{k)e 



Ei X b,B 



(2tt) 3 / 2 



e k@E ^H<,E,x\,B) 



A.n ) e idl ^ e * d3 ^y(xl, Ej xl B ) 



= V (| X l a E — X\ B I) -k a ,E ^( X a,E J X 1,b) ■ 

where V(k) is the Fourier transform of V(|x|) and 



(20) 



V(\x a ,E - x b,B\) *a,E ^( x l t Ei X \,b) 

= V(\x a , E - x b , B \) e d% ^ eijdJ ^*(xi„,xi„) 



(21) 
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is the Gronewold-Moyal or ^-product [T[JJI2n]- Notice that x l a E does not commutes with 
Pa, e but, however, it does commute with O^p^ E due to the antisymmetric character of 
0j§. This observation is of the outmost importance for arriving at Eq. (|2()jl . Furthermore, 
since only the electron coordinates are sensitive to the ^-product we drop, from now on, 
the subscript E in this particular symbol. We single out 



d 3 s0 1 (x)*0 2 (x) = J d 3 x0!(x)0 2 (x), (22a) 

J d 3 a;0i(x) *</> 2 (x) *0 3 (x) = J d 3 a; 3 (x) * ^(x) ★ 2 (x) 

d 3 x0 2 (x)*0 3 (x)*0 1 (x), (226) 



as the properties of the ^-product [H10IBlin] which will play a relevant role in our future 
developments. 

We now address to the problem of computing TCeb- By substituting Eq.(j20J) into 
Eq.()18jl one obtains, 



Heb*(x 



N 

c,b.^,b) = ^b,Bn(x b)B ) V 

a=l J 

X 4,b) 



6,B 



iV 

I 7 



^ A d 3 x 6)B v(|< )B - 4 >B |) * a *i 

a=l ^ 



d 3 zV(|;z|) 



*a ^(^c.E J X oI,b) 



a=l 

which in view of the arbitrariness of ^{x l c E , x l d B ) amounts to 
T~Leb = H E b > 



(23) 



(24) 



as an operator identity. Thus, the noncommutativity of the electron position observables 
does not affect the Hamiltonian Heb- This is a consequence of the continuous structure 
assumed for the ion background. 

It remains to study the modifications induced by the noncommutativity on He- 
One may convince oneself that 



'He = H Q + Ve , 
where Ho is given in Eq. (j8aj) . while 



(25) 



feiAj 

fciAi k2\'2 k^Xs 



k2^2 



k 1 X 1 k 2 X 2 



V 



Qa,E — Qb,E 



C fc 4 A4 C k 3 \ 3 



(26) 
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Through standard manipulations [13] , the right hand side of Eq. can be written 



&4A1 &2A2 



V 



Qa,E — Qb,E 



V- 2 [v\ 1 (aY®VxM] [Vx 3 (a)®VxM / d 3 x a , E / d^e"^— 



x e 



-i k2-x b 



V(| - f 6jB |) * a * b e 



where 



$ Sa (5) = (x|fcA) = V^e^A, 



(27) 
(28) 



is the free electron wave function, with 





1 







m = 





> m = 


1 



Furthermore, 



27m,- 



, i — 1, 2, 3 and = ±1, ±2, . . . 



(29) 



(30) 



is the periodically quantized momentum. The use of Eqs. (j22qJ) and (|225|) enable us to 
find 



k 1 X 1 k 2 X 2 



V 



Qa,E — Qb,E 



V 2 ^AiA 3 ^A 2 A 4 



k 3 X 3 /E4A4 

/d««..,/dv.v(|3:,,-4,.i) 



X 



= ^ 1 ^ 1+ fe,^/A 1 A 3 5A 2 A 4 e i ^ Afcl + ^) y d s sV(|3|)e^- fa \ (31) 

where the wedge product stands for 

kAp = k i Q%p>. (32) 
For arriving at the last term in the right hand side of Eq. (J31|) we took advantage of 

Q ix-k ^ e -ix-p = e ikAp e ig.k e -ig.p (33) 

The replacement of Eq. (|31|) into Eq. (J26)) yields 



2 °fciAi 
fciAi fc 2 A 2 £3-^3 &4A 4 

x / d 3 zV(|^)e iF - ( ^-^ 2) 



J V^ 1 5- — . A, , A, , i(fc 3 Afci + fc 4 Afc 2 ) 
fc 2 A 2 %+fc 2 , fc 3 +fc 4 ° A i A 3 °A 2 A 4 e 



fc 4 A 4 fc 3 A3 



2 V 5^ S 



d^VflzDe* 



-igA(fc-p) t J r r ^ 

h+q,Xl P-9, A 2 °P A 2 L fcAi 



(34) 
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This is the desired form of Ve in terms of creation and annihilation operators. It exhibits 
explicitly the noncommutativity. As is common practice, we have chosen q to designate 
the momentum transfer of the reaction p + k — > {p — q) + (k + q). 

As in the commutative case ^3], the contribution of the q = mode in the right 
hand side of Eq. (p£i|) cancels out those arising from Hb and Heb- This means that the 
noncommutativity does not destroy the electric neutrality. Hence, the whole modified 
system collapses into 

H = H + Hi , (35) 

where H is given by Eq. fJHaj) . while Hi reads 

2yj- g2 _^ _^ g— ig*A(fc — p) 

n i = — E ' E — ^ — 4 +3 -, Ai 4-?, a 2 % a 2 ^ Al . (36) 

At this point a digression is in order. Notice that, in contradistinction with the 
relativistic case, the commutative limit (Q E — * 0) in Eq. (j36|) exists and is well defined. 
To put it differently, the UV/IR mechanism [23], that contaminates noncommutative 
relativistic field theories, does not presently arise. This is of course due to the absence 
of ultraviolet divergences in the nonrelativistic case. It is a rather simple exercise to 
verify that Hi is Hermitean, as it must. 

We have so far developed the tools to compute some of the physical effects induced 
by the noncommutativity in the electron gas. We focus on the ground state energy 
eigenvalue and employ, as in the commutative situation, time-independent perturbation 
theory. We start by writing 

I (E<°> H, Ef> ) 2 

where {Ef*} are the excited states of H . Since H does not feel the presence of 
noncommutativity, its eigenstates and corresponding eigenvalues remain unchanged. 
Therefore, Eq.® still holds true. 

What come next is the computation of 8^ which, according to Eq.(|37|). reads 



£o = ^ 0) + (Eg 



Hi 



S, 



(i) 



E 



(0) 



"Hi I E Q 

2vre 2 



(0) 



V 



E'E 

kpq A l A 2 



-iq A (k — p) 



E 



(o) 



k+q 



fjAl C p-q,\2 C P. A 2 °k,\i 



E 



(0) 



• (38) 



As already indicated, the mode q = does not contribute to the right hand side of 
Eq. (|3*5j) . Then, straightforward manipulations lead us to 



E, 



(0) 



C l+q, Ai C P-€, ^2 C P' A 2 °k , X 1 



E 



(0) 



<(k F -p)Z(k F -k)5. 



p-q 



k ^AiA 2 



(39) 
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Observe now that, when replacing Eq.(|39|) into Eq.(|38|). the factor kills all 

noncommutative effects and, therefore, 



(i) 



E 



(i) 



(40) 



The computation of 



(2) 



(2) 



l(4 0) 




£f ) 


E 


(0) 





(41) 



is cumbersome. We shall not pause here to present the details but merely mention that 
it turns out to be given by 

-.2 

,(2)6 , 



(2) 



2a r 



N 



,(2)r 



+ e 



(6) 



(42) 



It is instructive to compare this result with its commutative counterpart, quoted in 
Eqs. flll|) and (|T2*|) . On the one hand, remains unaffected by the noncommutativity 
while, on the other hand, the exchange correlation energy term, e^ b , is modified as 
follows 



4 2)b -4 2)6 (©) 

= _3_ rdhi r d 3 k r d 3 f(i-k)e(i- P )e-"^^(i-a 

167T 5 7 q 2 7|k + q|>i 7|p-q|>i [q • (q + k - p)] (q + k - p) 2 

One may easily verify that e^ b (O) is real, as demanded by the hermiticity of Tii. This 
allows the replacement of the exponential by its real part, i.e., 



4 2)fe (6) 



3 /-d 3 q r , r , C(l-k)e(l-p) cos 2A; 2 qA(k-p) 

' ' d d k / d^P' 



167T 5 y q 2 Jk +fJ |>i y|p-q|>i [q-(q + k-p)](q + k-p) 2 



(44) 

Needless to say, the equivalence between Eqs.()43p and (jUj) can also be checked by 
direct computation. Since the argument of the trigonometric function depends on kp 
and, therefore, on V, the exchange correlation energy is no longer a constant. As 
consequence, the thermodynamics properties are modified by the noncommutativity. 

We have not been able of computing analytically the integral in Eq. (j44|) . To proceed 
further on, we assume that the global features of the system are insensitive to the 
direction of the vector Q % = l/2e^ k Q-' E k . We may, then, replace the right hand side of 
Eq. (pl4*|) by its average over all possible directions of 9. This is effectively achieved by 
integrating over the angles of 6 which yields 
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ef\e) s -iy dtv«r0 



^ d-k/" d*p «!- k )«!-P) 



16vr 5 7 q 2 7|k+q|>i 7|p-q|>i [q- (q + k - p)](q + k - p) 2 
x sin(fc|g|qx (k-p)|) 
fc 2 #|qx (k-p)| 



Let us concentrate on analyzing different limiting cases. For 9 = \9\ =0 one returns 
unambiguously to the commutative model. On the other hand, when 9 — > oo 6$ — > 
implying that £q < Eq . However, thermodynamic quantities, such as pressure and 

(2) (2) 

bulk modulus, will remain unaltered because the difference Sq — Eq is just a constant. 

The next step consists in bringing into play the Seiberg-Witten jl] map. By 
expanding the trigonometric function in Eq.(|45J) around 9 = one arrives at 

e^\6) = e^ b -J- 5 k*Re> + 0(9±), (46) 

where 

r=(*±[ d»k/ d 3 p ;(i-k);(i-P)igx(k-p)p (4?) 

J q i|k+q|>i J\p-q\>i [q • (q + k - p)](q + k - p) 2 

The convergence of the k and p integrals is secured by the fact that they run over finite 
intervals. On the other hand, power counting tell us that the improper q integral also 
converges. Hence, R exists and is well defined. The situation changes drastically for 
those integrals that act as coefficients of higher orders in 9. There, power counting 
indicates that they are divergent. The way out from the trouble consists in carrying 
out the q integral between and A being A a cutoff such that, as 9 — >• 0, l/k 2 F 9 goes to 
infinity faster than A. 

By collecting all the results, Eq.(pT7J) yields 



So = 

2a 



2.21 0.916 



En — N 



m 



+ 0.0622 In r. 
3 



k%e 2 R9 2 



0.094 



0(9 4 ) 



-^-k 4 F 9 2 R 

32tt 5 f 



0(9\r s \nr s 



(4£ 



h 32vr 5 

The noncommutativity certainly modifies the ground state energy and, as consequence, 
the ionization potential of the material being treated as an electron gas. Moreover, from 
the comparison of Eq. (J48|) with the commutative electron gas at non zero temperature 
[T7] one may conclude that 9 acts as a reference temperature [TO since d£ /d9 is a 
linear function of 9 much as the specific heat at constant volume is a linear function of 
the absolute temperature. This is the main outcome of the present paper, namely, the 
noncommutativity of the position observables lifts the degeneracy of the model and can 
be interpreted as if the electron gas would be at nonzero temperature. 
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